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^ ! Abstract 

C^ . A study of the i^-jet mass spectrum in e~^p -^ vX events at center-of-mass energy 

t^ , 300 GeV has been performed with the ZEUS detector at HERA using an integrated 

luminosity of 47.7pb~ . The mass spectrum is in good agreement with that expected 
from Standard Model processes over the i^-jet mass range studied. No significant excess 
attributable to the decay of a narrow resonance is observed. By using both e'^p -^ e^X 



•/^ ' and e^p -^ uX data, mass-dependent limits are set on the s-channel production of 

Qs ' scalar and vector resonant states. Couplings to first-generation quarks are considered 

O ■ and limits are presented as a function of the e^q and uq branching ratios. These limits 

f^ ■ are used to constrain the production of leptoquarks and i?-parity violating squarks. 
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1 Introduction 

A number of extensions of the Standard Model of elementary particles predict the ex- 
istence of electron-quark resonant states at high mass. Such states include leptoquarks 
(LQs) [|1| and i?-parity violating (/?p) squarks [^. The corresponding production processes 
could give a large cross section for high-mass i^-jet or e^-jet events. 

This paper presents an analysis of the ZEUS data aimed at searching for high-mass scalar 
and vector resonant states decaying into an antineutrino plus a jet. A similar search in 
the e"^-jet final states with the ZEUS data was published previously p. To avoid the 
constraints from a specific model, minimal assumptions are made about the properties of 
the resonant state. 

This analysis uses events whose observed final state has large missing transverse mo- 
mentum and at least one jet. These event characteristics correspond to an outgoing 
antineutrino and a scattered quark in e+p -^ VX scattering. The data-selection and 
event-reconstruction techniques are similar to those used for measuring the charged cur- 
rent (CC) cross section [^. The z/-jet invariant mass is calculated from the energies and 
angles of the final-state antineutrino and jet: 



Mj. =2E,E,et(l-cosO (1) 

where Ey and Ej^t are the energies of the scattered antineutrino and jet (assumed mass- 
less), respectively. The angle ^ is the laboratory-frame opening angle between the jet 
and the antineutrino. Since the antineutrino escapes detection, its momentum is deduced 
from all observed final-state particles by assuming conservation of energy-momentum in 
the event. 

In the following sections, the expectations of antineutrino-jet final states from the Stan- 
dard Model (SM) and from models that predict resonant states are first reviewed. After 
a summary of experimental conditions and data selection, the analysis is described and 
the reconstructed mass spectrum is presented. Since there is no evidence for a narrow 
resonance in either the z/-jet or the previously published e"'"-jet mass spectra, limits are set 
on the production of positron-quark resonant states using both data sets. The application 
of these limits to LQ and squark production is then discussed. 



2 Signal and background expectations 

High-mass z/-jet final states can be formed either through SM mechanisms or via processes 
that produce lepton-quark resonances. Figure |T| shows scattering mechanisms producing 






Figure 1: Processes with z/-jet final states in e~^p collisions. A scalar (S) or vector 
(V) intermediate state can be formed via a) s-channel or b) w-channel exchange. Weak 
charged current scattering c) forms the primary background to these processes. 



such final states in e~^p collisions. The CC scattering mechanism shown in Fig. |l|c forms the 
primary background in this search. Neutral current (NC) and photoproduction processes 
form negligible backgrounds since neither produces events with a large observed final-state 
momentum imbalance. 



2.1 Standard Model expectations 

The kinematic variables used to describe the process e~^p — >■ VX are: 

Q' = -l" = -(k - k'f (2) 

where P is the four-momentum of the incoming proton, and k and k' are the four-momenta 
of the incoming positron and the outgoing antineutrino, respectively. These variables 
are related by Q^ = sxy. The quantity x is interpreted as the fraction of the proton 
momentum carried by the struck quark, and y measures the fractional energy transferred 
by the W in the CC process. 

Assuming no QED or QCD radiation, the mass of the Vq system is related to x via 

M^ = sx (5) 

and the scattering angle, 6'*, of the outgoing antineutrino relative to the beam positron, 
as viewed in the Vq center-of-mass system, is related to y via 

cos e* = l-2y . (6) 

In leading-order electroweak theory, the CC cross section can be expressed as: 



(Pa^^{e+p) sGl ^ ^^^ 



dxdy 4:71 



(sfe) \y^^^-y-^^^^-y'^n m 



where Gp is the Fermi constant, Mw is the mass of the W boson, and Y± = 1±{1 — y)"^. 
The proton structure functions -F^*" and xFf^ , in leading-order (LO) QCD, measure 
respectively sums and differences of quark and antiquark parton momentum densities 0. 
The longitudinal structure function, -F^*^, contributes negligibly to this cross section 
except at y near 10. In the region of high mass (x — * 1) the structure functions F2'" 



and xF^^ are dominated by the valence quark distributions in the proton. For e+p 
colhsions, the scattering from down quarks dominates the cross section. The CC cross 
section peaks at small y, which leads to a cos 6'* distribution rising toward cos 6'* = 1. 

The largest uncertainty in the CC cross-section prediction arises from the parton densities 
of the proton. The parton density functions (PDF) are parameterizations which, at high x, 
are determined primarily from measurements made in fixed-target deep inelastic scattering 
(DIS) experiments. In the high-mass range (x ~ 0.6 corresponding to a z7-jet mass of 
230 GeV), the PDFs introduce an uncertainty of ~ 25% in the predicted e+p CC cross 
section 0. It should be noted that recent studies of PDFs suggest that the rf-quark 
density in the proton has been systematically underestimated for x > 0.3 [@, 0, |], 3- As 
an example, Yang and Bodek |^ propose a correction to the d/u quark density ratio in 



the MRS(R2) PDF ^■. 

5(^=Q.lx{x + l) (8) 

which fits the available data better. When this correction is applied to the CTEQ4D 



PDFs |T^, the increase in the predicted CC cross section (and the corresponding number 
of high-mass z7-jet events) ranges from 1.0% at x = 0.1 to 60% at x = 0.6. More recent 
PDF parametrizations P, ||, |12|, agree well with the corrected CTEQ4 for x up to 0.7. 



2.2 High- mass resonant states 

If a high-mass resonant state were produced at HERA, it could have a final-state signature 
similar to NC or CC DIS. Electron-quark states which couple to a single quark genera- 
tion and preserve lepton flavor are considered here. For e^p scattering, first-generation 
couplings of the form e^M, e^rf, e^u and e^d can be defined. 

These states are classified using the fermion number F = L + 3B, where L is the lepton 
number and B is the baryon number of the state. The coupling of positrons to quarks 
(e+M and e+c?) requires F = and the coupling of positrons to antiquarks {e~^u and e'^d) 
requires F = —2. In e~^p scattering, the F = states couple to the valence quarks of the 
proton and, for the same coupling, would have a significantly larger cross section than 
would the F = —2 states. 

Table |I] lists the 8 scalar and vector resonant states considered here, along with their 
charges and relevant decay modes. The e~^u and e~^d states would produce both e~^q and 
uq final states, which correspond to NC and CC event topologies, respectively. The other 
states would decay only to e'^q since a Vq mode would violate charge conservation. Some 
physics models incorporating high-mass resonances predict additional decay channels with 
final-state topologies different from DIS events. The branching ratios of each resonance 
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Table 1: Possible first-generation scalar and vector resonant states in e~^p scattering. The 
top half of the table lists color-triplet states with fermion number F = L + 3B = 0, 
while the bottom half lists those with F = —2. The left and right sets of columns list 
scalars and vectors, respectively. The e+d and e+u states can decay to both Pq and e+q. 
For the other states, only e+q decays are allowed since a Pq decay would violate charge 
conservation. 

into e+q, Pq and other final states are treated as free parameters except when specific 
models with restricted branching ratios are considered. 

In general, high-mass states formed by e+p collisions can have a combination of left- 
(Al) and right- (Xr) handed couplings. Because decays to right-handed antineutrinos 
must occur through left handed couplings, only left-handed coupled states {\r = 0) are 
considered for Pq decays. 

If a state with mass M^+g < i/i exists, the s-channel mechanism (Fig. la) would produce 
a resonance at Mj^j = M^+q in Pq decays. Additional contributions to the e+p cross section 
come from w-channel exchange (Fig. lb) and the interference with W exchange (Fig. Ic). 
The total e+p — > PX cross section with a resonance contribution can be written as pl 



d^a{e+p) _ d^a 
dxdy 



cc 



^2^Int 



CF. 



+ 



u/CC 



dxdy dxdy 



'^^^Ifcc d'^au d'^a. 



dxdy dxdy dxdy 



(9) 



The first term on the right-hand side of Eq. (|^) represents the charged current contribution 
from the SM. The second (third) term is the interference between the SM and w-channel 
(s-channel) exchange, and the fourth (fifth) term represents the w-channel (s-channel) 
exchange alone. The contribution of a single vector or scalar state has two free parameters: 
Me+g, the mass of the state and A, its coupling to e"'"-quark. The cos^* dependence of the 
state varies strongly for the different terms: it is uniform for a scalar state produced in the 
s-channel or a vector state produced in the w-channel, while it varies as (1 + cos 6'*)^ for 
a vector state produced in the s-channel or a scalar state produced in the u-channel M. 
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Table 2: First-generation leptoquark species considered in this analysis. The superscript 
L denotes chirality, while the subscript indicates the weak isospin. The electric charge, 
the production channel, and the allowed decay channels are also displayed. For positron 
beams, the charge changes sign, the helicity of the lepton is reversed, and the quarks and 
anti-quarks are interchanged. 

For the small couplings considered here, and if Mg+g < ^/s, the narrow resonance produced 
by the s-channel exchange would provide the dominant additional contribution over the 
SM background. The width of the s-channel resonance is given, e.g., for the Se+q with 
50 % branching to Vq, by 



M 



e^q 



e^q 



IGvr 



(2A2 



(10) 



so that if A^ is sufficiently small, the production cross section can be approximated by 
integrating over the s-channel contribution to the cross section. This leads to the narrow- 
width approximation for the total cross section of a single state |jl|: 



a 



NWA 



77 



4s 



2 ^ 



(11) 



where g(xo, M^+ ) is the initial-state quark (or antiquark) momentum density in the proton 



e~^q 

evaluated at Xq = M^+^/s and at a virtuality scale of M^+ , and J is the spin of the state. 



C'ql 



In the limit-setting procedure (Sect. ^, this cross section was corrected for expected 
QED and QCD radiative effects. The effect of QED radiation on the resonant-state cross 
section was calculated and was found to decrease the cross section by 5 — 25% as M^+q 
increases from 100 —>■ 290 GeV. For scalar resonant states, the QCD corrections |]13| 
raise the cross section by 20 — 30% for F = resonances. For F = 2 states, the QCD 
corrections lower the cross section by 5 — 30% in the 200-290 GeV mass range. No QCD 
corrections were applied to vector states because the calculation for such states is not 
renormalizable |T^ . 
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Table 3: Squarks predicted by SUSY that have flp decays into e^-jet or z/-jet final states. 
Listed are the squark production mechanism and decay channel. The k and j subscripts 
indicate the squark generation. Also shown is the corresponding resonant state from 
Table |1|. The decay modes with a Xi ' ^ire the i?-parity-conserving decay modes which 
produce neutralinos (x°) and charginos (xt)- These undergo further decays into SM 
particles. 

3 Resonant-state models 

In the absence of a clear resonance signal, limits can be placed on the production of states 
in models which predict a high-mass positron-quark resonance decaying to e~^q or Pq. 
Two such models are considered: (1) leptoquark (LQ) states with SU{3) x SU{2) x ^(1) 
invariant couplings and (2) squark states found in i?-parity violating supersymmetry 
(SUSY) models. 



3.1 Leptoquarks 

For SU{?)) X SU{2) x f/(l) invariant LQ couplings, there are 14 possible LQ species [Q. 
Such leptoquarks have no decay channels other than e^q or Vq. Table |^ lists those which 
have equal branching ratios into e+g and Vq decays. These scalar and vector LQ species 
correspond to the Se+n and Ve+d resonant states, respectively, with branching ratios fixed 
to /3e+g = /3p, = 1/2. 

3.2 SUSY 



In SUSY, conservation of baryon and lepton number is expressed in terms of i?-parity, 
Rp. It is defined as Rp = (— i)3-B+i+25^ where B is the baryon number, L is the lepton 
number and 5* is the spin of the particle. Ordinary SM particles have Rp = +1 while 
their hypothetical supersymmetric partners have Rp = —1. In versions of the theory in 
which i?-parity is not conserved, squarks (the SUSY counterparts to quarks) have the 




a) 



b) 




c) 



d) 



Figure 2: Lowest-order s-channel diagrams for first-generation squark production in 
e~^p collisions at HERA. Diagrams a) and c) are the /?p decays for dk and Uj squarks, 
respectively. The i?p-conserving decays are shown in b) and d). The decays of the 
charginos and neutralinos, x? and xt ^ i^^o SM particles depend on the parameters of the 
SUSY model and are not shown. 



same production mechanism as a generic scalar resonance. The squark flavors listed in 
Table || have flp decays into lepton-jet final states. Figures |a and c show the s-channel 
diagrams for these squark decays. The Uj and the c4 squarks behave like S^+d and S^+u 
resonant states, respectively (see Table ^, and the subscripts j and k denote the squark 
generation. Three generations are possible, but it is assumed that only a single generation 
has non-negligible coupling. These squarks would also be expected to have /2p-conserving 
decays into neutralinos (x°) and charginos ixt) (Figs. ||b and d) with multi-jet signatures 
different from e"'"-jet and z7-jet. A detailed discussion of these states, whose properties 
depend on many SUSY parameters, is beyond the scope of this paper. The branching 
ratios of squarks into e+-jet and P-jet, as well as other final states, are therefore treated 
as free parameters in this paper. 



4 Experimental conditions 

During 1994-97, HERA collided protons of energy Ep = 820 GeV with positrons of energy 
Ee = 27.5 GeV. The integrated luminosity of the data is 47.7pb^^. A detailed description 
of the ZEUS detector can be found elsewhere |jl5[. The primary components used in the 
present analysis are the central tracking detector (CTD) positioned in a 1.43 T solenoidal 
magnetic field, the uranium-scintillator sampling calorimeter (CAL) and the luminosity 
detector (LUMI). 



The CTD |jT6[ was used to establish an interaction vertex with a typical resolution of 3 cm 
in the beam direction for events considered in this analysis. Energy deposits in the CAL 
[ P^ ] were used to measure the positron energy and hadronic energy. The CAL has three 
sections: the forward^, barrel, and rear calorimeters (FCAL, BCAL, and RCAL). The 
FCAL and BCAL are segmented longitudinally into an electromagnetic section (EMC) 
and two hadronic sections (HACl, 2). The RCAL has one EMC and one HAC section. The 
cell structure is formed by scintillator tiles. The cells are arranged into towers consisting 
of 4 EMC cells, a HACl cell and a HAC2 cell (in FCAL and BCAL). The transverse 
dimensions of the towers in FCAL are 20 x 20 cm^. One tower is absent at the center 
of the FCAL and RCAL to allow space for passage of the beams. Cells provide timing 
measurements with resolution better than 1 ns for energy deposits above 4.5 GeV. Signal 
times are useful for rejecting background from non-ep sources and for determining the 
position of the interaction vertex if tracking information is unavailable. 



Under test beam conditions, the CAL has a resolution of 0.18/ y^E(GeV) for positrons 



^ The ZEUS coordinate system is right-handed with the Z axis pointing in the direction of the proton 
beam (forward) and the X axis pointing horizontally toward the center of HERA. The polar angle 6 
is defined with respect to the Z axis. 
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hitting the center of a calorimeter cell, and 0.35/ ^^/E(GeV) for single hadrons. The 
events of interest in this analysis have only hadronic jets, which impact primarily in 
the FCAL. In simulations, the jet energy resolution for the FCAL is found to average 
a/E = 0.55/v/^(GeV) © 0.02 §. 

To reconstruct the hadronic system, corrections were applied for inactive material in 
front of the calorimeter. The overall hadronic energy scales of the FCAL and BCAL are 
determined to within 2% by examining the Pt balance of NC DIS events [0. 

The luminosity was measured from the rate of the bremsstrahlung process e~^p — > e^p7 
19|] , and has an uncertainty of 1.6%. 



A three-level trigger similar to the one used in the charged current analysis was used to 
select events online [ffl. 



5 Event simulation 

Standard Model CC events were simulated using the HERACLES 4.6.2 EO] program 



with the DJANGO 6 version 2.4 [^] interface to the hadronization programs. First- and 



second-generation quarks are simulated, while third-generation quarks were ignored p2 
because of the large mass of the top quark and the small off-diagonal elements of the 
CKM matrix. The hadronic final state was simulated using the MEPS model in LEPTO 



6.5 [23], which includes order-a^ matrix elements and models of higher-order QCD radi- 



ation. The color-dipole model in ARIADNE 4.08 ||2^ provided a systematic check. The 



CTEQ4D parton distribution set ||TT| with the Yang-Bodek correction, Eq. (g|), was used 
to evaluate the nominal CC cross section, and the unmodified CTEQ4D PDF was used 
as an alternative PDF with smaller d-quark density. 



Simulated resonant-state events were generated using PYTHIA 6.1 ||2^. States with 
masses between 150 and 280 GeV were simulated in 10 GeV steps. This program takes into 
account the finite width of the resonant-state, but only includes the s-channel diagram. 
Initial- and final-state QCD radiation from the quark and the effect of LQ hadronization 
before decay are taken into account, as is initial-state QED radiation from the positron. 

Generated events were input into a GEANT 3.13-based simulation [^] of the ZEUS 



detector. Trigger and offline processing requirements as used for the data were applied 
to the simulated events. 
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6 Event selection 

Events were selected with cuts similar to those used in the CC cross-section measurement 
from the same data ||^. The events were classified first according to 70, the hadronic 
scattering angle of the system relative to the nominal interaction point [Q. If 70 was 
sufficiently large, i.e. in the central region, tracks in the CTD were used to reconstruct 
the event vertex. On the other hand, if 70 was small, i.e. in the forward region, the 
hadronic final state of such i^-jet events was often outside the acceptance of the CTD, 
and thus the vertex position was obtained from the arrival time of particles entering the 
FCAL. The following selection cuts were then applied: 

• to select high-mass PX states, events were required to have substantial missing trans- 
verse momentum: f^ > 20 GeV; 

• a cut of y < 0.9 discarded events in which the kinematic variables were poorly recon- 
structed; 

• events with ^t/Et < 0.4 (where E^ denotes the total transverse energy measured 
in the event) were removed to reject photoproduction background. For events with 
7o < 0.4, this cut was increased to 0.6; 

• NC background was removed by discarding events with identified positrons; 

• non-ep collision events caused by beam-gas, halo muons, and cosmic rays were removed 
by a series of standard cuts based on the general topology expected for events from ep 
collisions originating from the interaction region at the correct beam-crossing time. 

The final sample contains 829 events. 

The momentum carried by the antineutrino is extracted from the /V and the longitudinal 
momentum variable {E — Pz) of the event; distributions are shown in Fig. ^. The data 
and SM predictions agree except for /V > 90 GeV, where a slight excess is observed in the 
data. The (E — Pz) distribution peaks near 10 GeV. These distributions are very different 
from those of NC events, which have small /V and an [E — Pz) distribution peaked near 
twice the positron beam energy. These differences arise from the undetected final-state 
antineutrino in this sample. 



Jets were identified using the longitudinally-invariant fcr-clustering algorithm |^ in in- 
clusive mode [^. At least one jet was required with transverse momentum P^ > 10 GeV. 
Fig. 1^ shows the distributions of the pseudorapidity, rj, of the highest P-f jet f\ Also shown. 



The pseudorapidity is defined as ry = — ln(tan(6'/2)). 
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Figure 3: (a) The fx distribution for the final event sample, (b) The {E — Pz) distribu- 
tion. The points represent the data and the histogram is the SM MC prediction. 
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Figure 4: Comparison of jet distributions in the data and Monte Carlo, (a) The pseudo- 
rapidity, rj, of the highest P^ jet in each event, (b) The energy of the highest P^ jet. (c) 
The missing transverse momentum /V in each event when the highest P^ jet is excluded. 
The points represent the data and the histogram the SM Monte Carlo prediction. 
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Figure 5: (a) The distribution of the energy of the final-state antineutrino in the lab 
frame, Eiy. (b) The distribution of cos^,y, where 6y is the polar angle of the scattered 
antineutrino in the lab frame. The forward direction (cos^j, = 1) corresponds to the 
proton beam direction. The points are the data and the line is the SM Monte Carlo 
prediction. 



for each event, are the energy of the highest P^ jet and the /V when the momentum of 
the highest P-f jet is excluded. Reasonable agreement is observed between the data and 
SM predictions in each case. 

The outer boundary of the inner ring of FCAL towers was used to define a fiducial cut 
for the jet reconstruction. The centroid of the jet with the highest P^ was required to 
be outside a 60 x 60 cm^ box on the face of the FCAL centered on the beam pipe. This 
restricts the pseudorapidity of the jet to be less than roughly 2.6. This requirement 
removes 25 events, bringing the total sample to 804 events. 
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7 Mass and 6* reconstruction 

It was assumed for the resonance search that all the missing momentum is carried away 
by one antineutrino. The invariant mass of the V-]ei system, M^j, was calculated using 
Eq. (P using only the highest Pf jet. The jet direction was determined from the vector 
formed by the event vertex and the jet centroid in the calorimeter. The neutrino energy 
and angle were calculated as: 



E, 



V 



COS 6„ 



Ft' 


+ {E 


-Pz)l 


2- 

Ft' 


{E- 
-{E 


Pz)u 

-Pz)l 



fr' + {E- Pz)l 



where {E — Pz)v = 2Ee — {E — Pz). Distributions of the reconstructed antineutrino energy 
and polar angle in the laboratory frame [E^, and cos 6,^) are shown in Fig. |^. Reasonable 
agreement is observed between data and the SM prediction. Monte Carlo simulations of 
resonant states indicate that the antineutrino energy and polar angle were measured with 
average resolutions of 16% and 11%, respectively. The average systematic shift in E^ was 
found to be less than 2%, while the shift in 6i, was less than 1%. 

Monte Carlo simulations of resonant states were used to determine the resolution and 
estimate the possible bias for the reconstructed mass. The mass resolution was obtained by 
performing a Gaussian fit to the peak of the reconstructed mass spectrum. For resonant- 
state masses from 170 GeV to 270 GeV, the average mass resolution was found to be 7%. 
The peak position of the Gaussian differed from the generated mass by less than 2% over 
the entire range. 

Note that energy-momentum conservation, assumed in order to calculate Ei, and 9,^, does 
not apply when undetected initial-state radiation (ISR) from the beam positron occurs. 
At high masses, QED radiation results in an underestimate of E^, and an overestimate 
of Oy. This, as well as final-state QCD radiation, results in lower reconstructed masses, 
leading to an asymmetry in the expected mass distribution. In a simulation of a resonance 
of mass 220 GeV, only 1% of events had an M^^j more than 20% higher than the true 
mass, while 16% had an M,yj more than 20% lower than the true mass. 

In contrast to the resonance search, setting cross-section limits on e~^p -^ vX processes 
requires that a specific production mechanism be assumed. For this reason, an invariant 
mass, M^js, was calculated using all of the jets in the event with P^ > 10 GeV and r/ < 3. 
Monte Carlo studies show that, for narrow resonant states, using multiple jets gives more 
accurate mass reconstruction for events with more than one jet (for masses above 150 
GeV, 12% of the simulated LQ events have multiple jets). 

The selection cuts described in Sect. ^ determine the kinematic region where mass re- 
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Figure 6: Acceptance in the cos 9*-M^j plane. The shaded areas are the regions excluded 
by the requirements of /V > 20 GeV, y < 0.9 and jet- containment assuming an eq — > Vq 
scattering at the nominal interaction point. No detector simulation is included. The 
dotted 7o = 0.4 line shows the boundary between events using the FCAL timing vertex 
(above) and the CTD tracking vertex (below). 
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Figure 7: The distribution of the final event sample in the Mj^j-cos 6* plane. Solid points 
indicate events reconstructed with a tracking vertex; open circles events reconstructed 
with a timing vertex. 



construction is possible. Figure ^ shows the approximate regions in the cos 6*-M^j plane 
which are excluded by the requirements of fx > 20 GeV, y < 0.9 and the jet contain- 
ment for events originating from the nominal interaction point. In the unshaded regions, 
acceptance is typically ~ 80%. The variable 7 denotes the scattering angle of the struck 
quark. Events above the 70 = 0.4 line typically use the FCAL timing vertex, while those 
below this line use the vertex found from CTD tracking. 
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8 Mass and cos^* distributions 

Figure ^ shows the distribution of events in the M^j -cos 9* plane. The events populate 
the region of large acceptance described in Fig. p. 



8.1 Systematic uncertainties 

The systematic uncertainties in the predicted rate of events range from about 7% at 
M^j ^ 100 GeV to about 20% at M^j ^ 220 GeV, and over 40% at M^j ^ 260 GeV. The 
major sources of these are uncertainties in the calorimeter energy scale, uncertainties in the 
simulation of the hadronic energy flow (established by comparing results from the nominal 
LEPTO MEPS model with a Monte Carlo sample using the alternative ARIADNE model) 
and uncertainties in the parton distribution functions. 

Potential sources of systematic error which were found to have negligible effects in- 
clude reasonable variations of the selection cuts, background-contamination uncertainties, 
timing-vertex uncertainties, and the uncertainty in the luminosity determination. 



8.2 Comparison with Standard Model 

In Fig. |^(a), the observed mass distribution is compared to the SM predictions from 
Monte Carlo simulations using the CTEQ4D parton densities [0 and the CTEQ4D PDF 
modified by the Yang-Bodek correction of Eq. (||). The predictions using the CTEQ5 P] 
or the NLO QCD fit by Botje are similar to the modified CTEQ4D predictions. For 
Myj > 180 GeV, the data tend to lie above the expectations. There are 30 events observed 
in this region, while 21.5 ± 3.3 are predicted (16.0 ± 2.4 events for CTEQ4D without the 
correction of Eq. (|^)). The uncertainty on the predicted number of events is due to the 
effects described above. 



Figure ^ shows the cos 9* distribution of the events with M^j > 180 GeV together with 
the distribution expected for decay of a narrow scalar resonance (normalized to 9 events). 
In the cos6'* < 0.4 region where the DIS background is suppressed, 8 data events are 
observed while 3.6 ± 0.5 SM events are expected. 

Given the limited statistics in the present data and the systematic uncertainties of the 
SM predictions, the observed mass spectrum is compatible with SM expectations. 
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Figure 8: (a) The mass distribution for the data (points) and Monte Carlo (histograms). 
The dashed hne shows the predicted mass spectrum when the CTEQ4D PDFs are used, 
while the solid curve shows the distribution predicted when the o?-quark density is en- 
hanced using the Yang-Bodek correction (see Eq. (D). (b) The ratio of the number of 
events observed to the number expected, N"''^ /N'^^^, obtained using the Yang-Bodek cor- 
rection. The shaded band indicates the systematic error in the SM expectation. The 
dashed line shows the SM expectation when the Yang-Bodek correction is not imple- 
mented. The error bars on the data points are calculated from the square root of the 
number of events in the bin. 
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Figure 9: The cos 6** distribution of events with M^j > 180 GeV. The dashed hne shows 
the predicted cos^* spectrum when the CTEQ4D PDFs are used, while the sohd curve 
shows the distribution predicted when the d-quark density is enhanced using the Yang- 
Bodek correction (Eq. (D). Also shown is the cos 6'* distribution for a scalar resonance 
(dotted line) normalized to 9 events. 
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Figure 10: a) The reconstructed mass spectrum using multiple jets for data (points) 
and SM expectation (histogram), b) The mass spectrum using multiple jets after the cut 
(cos 9* < cos 6'j^ax) fo^ the scalar resonance search has been applied. 
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9 Limits on resonant-state production 

Since there is no evidence for a narrow resonance in the i^-jet data, hmits may be set on 
the production of the resonant states hsted in Table |I[ Since such states would need to 
have a positron as well as an antineutrino decay channel, the cross- section limits were set 
using these P-iet data along with the e"'"-jet data previously reported |^. Only couplings 
A < 1 are considered. The limit-setting procedure assumes the states have the same 
production and decay mechanism as the Monte Carlo used to generate the resonance 
events. The invariant mass reconstructed using the neutrino and all jets with /V > 10 GeV 
and rj < 3.0, Mj^j^, was used to set limits. The mass spectrum reconstructed with this 
technique is shown in Fig. |lO|a, and is similar to that from single jets (Fig. H). 



The limit-setting procedure requires two parameters at each value of M^jg-. the mass 
window, AMujs, and an upper cut {cos 9^^^) on the measured value of cos^*. Simulations 
of both SM background and resonant signals were used to find values for these parameters 
which optimize observation of a signal relative to DIS background. For a scalar resonance 
with a V-iet final state, AM^jg ranged from 20 to 35 GeV in the 160-280 GeV mass range, 
while in the same range cos^j^^^^ increased from 0.2 to 0.8. For a vector resonance in the 
same M^jg range, AM^jg increased from 15 to 35 GeV, while cos 9^^^ increased from 0.6 
to 0.84. The mass spectrum after applying the optimal cos^* cut for the scalar search is 



shown in Fig. |TO|b. A similar optimization procedure, performed for the e+-jet final state 



using the NC data, has been described in a previous publication 0. 

To find the 95% confidence level (CL) upper limit on the resonant-state cross section, 
(Jiim, a likelihood is calculated using the Poisson probability for each decay channel: 



where C is the luminosity, (3c is the branching ratio of the decay channel, N°^^ is the 
number of observed events, N^^^ is the expected number of DIS background, and ec is the 
acceptance calculated from resonance Monte Carlo. The subscript c denotes the decay 
channel, which for this analysis is either uq or e^q, for the CC-like and NC-like final states, 
respectively. If more than one channel was used to set a limit, the likelihoods for each 
channel were multiplied together to get the total likelihood, L{a). A flat prior probability 
density for the cross section a was assumed, such that the probability density, /(cr), is 
simply /(cr) oc L{a). A limit was then obtained on the cross section, crum, by solving: 

/"CTlim rOO 

/ daf{a) = 0.95 / daf{a) (13) 

Jo Jo 
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and the resulting cross-section limit was converted to a coupling limit Aum using the NWA 
(Eq. (|llD). Note that using two channels does not always produce a stronger limit than 
using a single channel. 

The limits on A depend on the accuracy of the NWA. Comparisons between the NWA 
and the full resonant-state cross sections show that the NWA was too high by up to a 
factor 1.7 for Sf.+u. This was corrected for in setting the limits. For all other states, the 
NWA provides a reasonable approximation of the full resonant-state cross section in the 
mass and coupling ranges studied. 



Figure |TT] shows the limits obtained for the four scalar resonant states of Table |T] as a 
function oi (3e+q and /5pg, the branching ratios into e^q and z/g, respectively. The equivalent 



plots for vector resonant states are shown in Fig. O. The limits were calculated for 



coupling strengths of A = 0.05 and A = 0.10, as well as for coupling A = 0.31 ~ v^lvra. 
For the e'^u and e+J resonances (a and d in Figs. |lT] and |12D, uq decays are forbidden by 
charge conservation, so the limits are set using only the e^q channel. The e^u and e^d 
resonances (b and c) can provide both e+g and Vq decays, so limits are calculated using 
the e^-jet and z7-jet data sets separately and combined. The combined e^q + Vq limits, 
which assume jS^q + P^+q = 1; are largely independent of branching ratio. The hmits 
obtained using only the e^-jet (or the z/-jet) data allow for decay modes other than e~^q 
and Pg, so the e~^q and the Pq limits are applicable to a wider range of physics models than 
the combined e'^q+Vq results. The systematic uncertainties on the predicted background 
described in Section 8.1 were found to change the excluded mass limits by less than 1% 
for M^j > 220 GeV, and have therefore been neglected. 

The e^q and Vq data have also been used to set limits on scalar and vector resonances 
with second generation quarks. Assuming a coupling strength of A = 0.31 the mass limits 
for e~^s states decaying with 50 % branching ratio to e^q and with 50 % to Pq are 207 GeV 
for a scalar and 211 GeV for a vector state. 



For comparison, the limits on scalar resonances obtained by the DO experiment [^ at the 
Tevatron are shown by the shaded region. These limits are independent of both coupling 
and quark flavor. Similar results to those presented here have been published by the HI 
experiment [^ . 



10 Model-dependent limits 

The limits on generic resonant states were converted to limits on the production of LQ 
and squarks that have e~^q and Pq decays. Figure |T3| shows the limit on the production 
cross section, a^^, for scalar and vector resonant states. Limits derived from e~^q [Pq) 
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Figure 11: The branching ratios into e^q and Pq (shown on the left and right axes, 
respectively) vs. excluded mass for the scalar resonant states listed in Table |l]. For each 
limit curve, the area to the left of the curve is the excluded region. Results for e~^u, e~^d, 
e^u and e~^d resonant states are shown for coupling strengths of A = 0.05, A = 0.10, 
and A = 0.31. The shaded region in each plot shows the mass range excluded by the DO 
experiment. For (a) e'^u and (d) e^d resonant states, limits were set using only e^q data 
since Vq decays are forbidden by charge conservation. The (b) e~^d and (c) e^u states have 
both e^q and Vq decay channels. The dotted line corresponds to only Vq data, the shaded 
line corresponds to only e'^q data, and the solid black line corresponds to both the e~^q 
and the Vq data sets. The combined limits were calculated assuming that (3pq + (3e+q = 1- 
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Figure 12: The branching ratios into e^q and Vq (shown on the left and right axes, 
respectively) vs. excluded mass for the vector resonant states listed in Table |l|. For each 
limit curve, the area to the left of the curve is the excluded region. Results for e^u, e~^d, 
e^u and e'^d resonant states are shown for coupling strengths of A = 0.05, A = 0.10, and 
A = 0.31. For (a) e^u and (d) e~^d resonant states, limits were set using only e^q data 
since Vq decays are forbidden by charge conservation. The (b) e~^d and (c) e~^u states have 
both e+g and Pq decay channels. The dotted line corresponds to only uq data, the shaded 
line corresponds to only e^q data, and the solid black line corresponds to both the e~^q 
and the Vq data sets. The combined limits were calculated assuming that /5pq + Pe+q = 1- 
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Figure 13: (a) Limits on the total production cross section for a narrow scalar resonant 
state, (b) The corresponding limits for a narrow vector resonant state. Limits derived 
from e~^q (Vq) assume a branching ratio Pe+q iPuq) of 1/2, while the combined e~^q + Vq 
limits assume branching ratios of P^+q = Puq = 1/2. 
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Figure 14: (a) The limits on the couphng Aum for an S",^ LQ. (b) The same for a V^ 
LQ. Results from the Vq and e^q channels are shown, along with the limits obtained by 
combining the two channels. Also shown is the limit for second generation LQ's (dashed- 
dotted line). In both plots, the horizontal line indicates the coupling A = 0.31 ~ v47ra. 
For comparison, representative limits from the Tevatron [^ and LEP [^ are also shown. 
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assume a branching ratio P^+q {Pvq) of 1/2, while the combined e^q + uq hmits assume 
branching ratios of /?e+g = Ppq = 1/2- 



10.1 Leptoquarks limits 

The cross-section hmits were converted to hmits on leptoquark couphng using Eq. (|ll|). 
Figure |l^ shows the couphng hmits for the Sq and Vg^ LQ species hsted in Table ||. If a 
coupling strength A = 0.31 ~ v47ra is assumed, the production of an Sq LQ is excluded 
up to a mass of 204 GeV with 95 % CL, while the production of a V^^ LQ is excluded 
up to a mass of 265 GeV. When the Vq and e~^q limits are combined, the resulting limits 



exclude approximately the same mass range as the e"'"g-only limit. Also shown in Fig. |n 
is the limit curve for second generation LQ's of the type V^ produced as an e~^s resonance. 
The combined limits from e~^q and Vq decays are shown. For comparison, limits from the 



DO experiment with a branching ratio of f3e+q = 1/2 are shown ||29[. Also included are 
LQ limits from the OPAL experiment at LEP 



10.2 SUSY limits 

Limits were set on the production of the squarks listed in Table ||. In addition to ^p decays 
into e~^q and Vq, squarks can also have i?p- conserving decays into other final states. To 
remove the dependence on the branching ratios into these i?p- conserving states, limits 
were set on the quantity X^/P, where /3 = Pe+q + Puq- The limit-setting procedure does 
not account for possible contributions to the e+-jet and I^-jet channels from i?p-conserving 
decays. Limits on dk and Uj are shown in Fig. ^. Because Pe+q = Pvq for the dk decays, the 
combined e^qWq limits are shown along with the limits obtained from the individual decay 
channels. For the Uj squark, /3 = Pe+q since Vq decays would violate gauge invariance. 
Previous limits on ^p-squark production from smaller data sets have been set by the HI 
experiment [^ . 



11 Conclusion 

A study of the z/-jet mass spectrum in e^p — > VX events at center-of-mass energy 300 GeV 
has been performed with the ZEUS detector at HERA using an integrated luminosity of 
47.7pb~^. Events with topologies similar to high-Q^ charged current DIS were selected. 
The invariant mass, M^j, was calculated from the jet with the highest transverse energy 
and the antineutrino four-momenta. The jet momentum was measured directly, while 
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Figure 15: The limits on the couphng X^/J3, where /3 = Pe+g + Puq- For the d^ squark, 
Pe+q = Puq, SO results from the P-jet and e"'"-jet channels are shown, along with the limit 
obtained by combining the two channels. The Uj limits are calculated using only e~^q data 
since the gauge invariance of the SUSY superpotential requires Ppq = 0. 
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the antineutrino inomentum was deduced from the energy-momentum imbalance mea- 
sured in the detector. No evidence for a narrow resonance was observed. This analysis 
complements an earlier search for narrow resonances in the e^-jet final state. 

In the absence of evidence for a high-mass resonant state, the e"'"-jet and z7-jet data sets 
were used to set limits on the production cross section of scalar and vector states decaying 
by either mode. Sensitivity to a resonant signal was optimized by restricting the center- 
of-mass decay angle to remove most DIS background and by choosing an appropriate 
mass window. The resulting cross-section limits were converted to coupling limits on 
e~^u, e'^d, e^u and e'^d resonant states. 

First-generation couplings between initial- and final-state quarks and leptons which con- 
serve flavor and electric charge were considered. Limits were calculated as a function of 
the e~^q and Vq branching ratios for small couplings and do not depend on a specific pro- 
duction mechanism. For resonances with both e~^q and Vq decays, using both the e"'"-jet 
and i^-jet data gave limits which are largely independent of the branching ratio if the state 
is assumed to have no additional decay modes. 

The limits on generic resonant states were used to constrain the production of leptoquarks 
and i?p-violating squarks. For leptoquark flavors whose branching ratios into e^q and Vq 
are the same, exclusion limits of 204 GeV for scalars and 265 GeV for vectors were obtained 
if a coupling strength A = 0.31 is assumed. Limits on the production of Uj and dk squarks 
were obtained directly from the limits on e^d and e~^u resonances, respectively. 
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